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Using non-linear system theory and numerical simulations we map out the static and dynamic 
phase diagram in zero applied field of a spin torque oscillator with a tilted polarizer (TP-STO). 
We find that for sufficiently large currents, even very small tilt angles (0 > 1°) will lead to steady 
free layer precession in zero field. Within a rather large range of tilt angles, 1° < f3 < 19°, we find 
coexisting static states and hysteretic switching between these using only current. In a more narrow 
window (1° < (3 < 5°) one of the static states turns into a limit cycle (precession). The coexistence 
of static and dynamic states in zero magnetic field is unique to the tilted polarizer and leads to large 
hysteresis in the upper and lower threshold currents for TP-STO operation. 



Spin torque, or the transfer of angular momentum from 
spin polarized electrons to magnetic moments l2[ , cur- 
rently receives an increasing interest due to potential use 
in magnetoresistive memory (MRAM) and in microwave 
signal_generators, so-called spin torque oscillators (STO). 
Hi, H 0, 0, @, 0, 0- While the first spin torque devices 
were based on (pseudo-)spin valves with in-plane mag- 
netizations, recent devices utilize perpendicularly mag- 
netized layers to achieve both higher stability in MRAM 
and zero-field operation in STOs. The resulting static 
and dynamic phase diagrams have been studied in detail 



In this Letter, we study the static and dynamic phase 
diagram of a spin valve where the magnetization of the 
fixed layer is tilted at an arbitrary angle out of the film 
plane. This so-called Tilted Polarizer STO (TP-STO) 
has the significant advantage of zero-field operation while 
maintaining a high microwave output signal [13]. We here 
show, using non-linear system analysis, that the deviation 
from in-plane orientation creates a surprisingly rich phase 
diagram with coexistence of different static and dynamic 
states within a certain range of the polarizer tilt angle /?. 
We determine the stable precessional states using mag- 
netodynamic macrospin simulations, and study the hys- 
teretic switching between both different static states and 
between static and dynamic states. The coexistence of 
static and dynamic states in zero field is unique to the 
TP-STO and disappears for polarizer angles outside of 
1° < (3 < 19°. As a consequence, the TP-STO can ex- 
hibit unexpected large current-driven hysteresis in both 
the upper and lower threshold currents for precession. 

Inset I in Fig. [T] shows the schematic structure of the 
TP-STO. While all films are deposited in the x-y plane, 
the fixed layer magnetization, M, lies in the x-z plane, 
with an angle (3 w.r.t. the x-axis. The time-evolution 
of the unit vector of the free layer magnetization rh 
is found from the Landau-Lifshitz-Gilbert-Slonczewski 
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FIG. 1: Phase diagram of static and dynamic states of a 
TP-STO with nodes (N), spirals (S) and Limit cycles (L). 
Parentheses denote states only found through dynamical sim- 
ulation. Inset I: The TP-STO structure: m and M are the 
magnetization vectors of free and fixed layer respectively. M 
lies in the x-z plane with angle /3 w.r.t. the x-axis. 6 and (j) are 
the polar and azimuthal angles of the free layer magnetization 
m. Inset II: Phase diagram of the high- J region. 



(LLGS) equation 0,0, 

-\j\rh x H e ff + OLrh x 



dm 
~dt 



dm 
~dt 



\j\ajm x (m x M), 

(i) 

where 7 is the gyromagnetic ratio, a is the Gilbert damp- 
ing parameter, and aj = cJ, with J being the electric 
current density and c a constant including material pa- 
rameters and fundamental constants. The electric current 
is defined as positive when it flows from the fixed to the 
free layer and normalized by Jo = 10 8 A/cm 2 . The effec- 
tive field tl e ff carries the contribution of an anisotropy 
(easy axis) field along the x-axis and a demagnetiza- 
tion (easy plane) field H^. It should be noted that the 
applied magnetic field is set to zero throughout this Let- 
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ter. For the results presented here, |7| = 1.9 • 10 11 Hz/T, 
a = 1CT 2 , HoH d = 1 T, [ioH k = 10" 2 T, with /x being 
the vacuum magnetic permeability. We use a symmetric 
torque term and a sinusoidal angular GMR dependence 
to not introduce further complications which would ob- 
scure the main results. 

Eq. (pQ) can be transformed into the following set of 
differential equations in spherical coordinate system: 
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(aU + W), 
aU-W 



sine 



(2) 



with U = [Hk cos 6 sin 6 cos 2 (j) + Hd cos 6 sin — 
aj cos /? sin (/>] , and W = [ — Hk cos sin sin <fi — 
aj (cos cos cos /? — sin sin /?) ] . 

By setting (U, W) = 0, we get a series of possible 
equilibrium solutions = §i(ft,J),^i = (j>i(/3, J), where 
i < i t is the ith solution of total i t solutions. However, 
these i t equilibrium states are not all stable. We linearize 
Eq.([2]) in the vicinity of (0$, fa) and get: 



'A(P,J,di,fc) B(P,J,di,fc) 



(3) 



where A, B, C and D are the explicit functions of vari- 
ables/3, J and other material parameters. Following 
Ref. [14], the eigenvalues of the corresponding Jacobian, 
which determine the stability of the system, can there- 
fore be solved and can always be expressed as: /ii 5 2 = 
E(ft, J) ± yj F\ft, J). For a solution to be stable, it must 
satisfy 5ft {/ii ,2} < 0. For real eigenvalues and F > 0, 
the eigenvalue with larger magnitude dominates and de- 
fines the only eigenvector governing the approach towards 
the final state, in this case a node (N). For F = the 
two eigenvectors are identical and again define a node. 
For F < the complex conjugate eigenvalues define two 
complex eigenvectors generating an oscillatory trajectory 
towards equilibrium, characteristic of a spiral-like (S) so- 
lution [in. 

As an illustrating example, we show the procedure for 
determining the stability and type of solution for the 
case ft = 0°, i.e. for a conventional in-plane spin torque 
MRAM cell in zero field. The well known solution for neg- 
ative current is (0, <p) = (tt/2, 0), i.e. parallel alignment of 
the free and fixed layer magnetizations. Expanding Eq. 
(PP) around this point yields: 



aj - a (H d + H k ) -H k - aaj 
H d + H k + aaj aH k + a j 



, (4) 



with eigenvalues: 
aH d 

Mi, 2 = — 



' oiH d 



H k (H k + Hd)(ai-l) + f(aj), (5) 



where f(aj) — aaj (Hd — 2H k — aaj). By entering the 
parameters into Eq. ([5]), we see that the type of solution 
and its stability depend upon the value of aj (i.e. J), 
if all other parameters are fixed. For \J/Jq\ < 1.5 the 
solution is of spiral type (£), while outside this region, 
where the torque is larger, the solution is a node (N). 
Point A in Fig. [T] denotes the S — > N transition at neg- 
ative currents. This result is well known in conventional 
spin torque switching where switching between S states 
proceeds by slow spiraling out of the unstable state and 
into the stable S state, while large currents will switch 
the magnetization without much precession into a stable 
N state Q. 

Following this procedure, we now construct the static 
part of the phase diagram in Fig. [T] by finding all eigen- 
values in the parameter space 0° < ft < 90° and \J/Jq\ < 
1.6 • 10 3 . While the entire parameter range was studied, 
Fig. CD focuses on 0° < ft < 30° and \J/J \ < 10, where 
coexistence of several different stable solutions are ob- 
served. It should be noted that only the static solutions 
(S and N) can be found from the eigenvalue analysis. 
To find the precessional states (L) we have to resort to 
numerical simulations below. However, according to the 
Poincare-Bendixson theorem [H, [l7| , the only possible 
final states of the system are either static states (fixed 
points) or limit-cycles (self-oscillation) and chaos is pre- 
cluded since the free layer evolves on the unit sphere sur- 
face (m = 1) 18, In regions where there are neither S 
nor TV solutions we can hence infer that steady precession 
(L) must take place. 

For small enough J there is a single S state at all tilt 
angles ft, corresponding to the usual P/AP orientation 
of the free layer with respect to the in-plane projection 
of the fixed layer magnetization. For ft > 2° this static 
state disappears with increasing current, and as discussed 
above, we can infer a precessional L state in this region. 
The crossover between S and L regions defines the crit- 
ical current for the onset of precession (J c i), where the 
negative damping from the spin polarized current desta- 
bilizes the S state and sustains continuous precession. It 
is hence possible to have zero-field TP-STO operation 
down to very small tilt angles if only large enough cur- 
rent densities can be realized. At yet smaller tilt angles 
(as in the in-plane case above) the eigenvalue analysis 
indicates that the S state changes into a node (N) with 
increasing current. It will however become apparent in 
the magnetodynamic simulations below that for ft > 1° 
this region also contains a precessional state (hence the 
additional label (L)). 

If the current is increased further, precession stops at 
an upper threshold current (J C 2) where the L state turns 
into a single spiral state located close to the north/south 
pole of the unit sphere. If we increase the current in the 
N/ (L) region the same L — > S transition occurs but in 
addition, the node also remains stable. We hence observe 
a region S/N where two different stable static states co- 
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current swept from (a) negative to positive and (b) positive 
to negative. Inset: the onset critical current density for pre- 
cession vs. ft. 

exist. The two states are located far from each other at 
two different points on the unit sphere: the north/south 
pole (5), and at P/AP alignment (N) respectively. As 
will become clear in the magnetodynamic simulation be- 
low, this separation allows both states to be realized by 
only sweeping the current. 

Finally, at large currents (inset II of Fig. [Tj) there is a 
node state at P/AP alignment for ft < 40° and a spiral 
state close to the north/south poles for ft > 40°. At ex- 
tremely high currents, the S state gradually turns away 
from the poles, approaches P/AP orientation, and finally 
replaces the N state at about J/ Jo > 10 3 as the free layer 
magnetization aligns completely with the fixed layer. At 
these current densities (J ~ 10 11 A/cm 2 ) any real sample 
would break down; we include this region for complete- 
ness. 

To determine the dynamic states and also study the 
hysteretic switching between the coexisting static states, 
we now solve Eq. (pQ) using numerical simulation within 
the macro-spin approximation. To simulate actual hys- 
teresis loops as a function of current, we start out at 
very large negative current, let the simulation reach a 
steady state, determine the type of state and its dynamic 
or static properties, and then let this state be the ini- 
tial condition for the next simulation at the next current 
step. At very large negative (positive) currents, the free 
layer always aligns (anti-aligns) with the fixed layer, as 
was confirmed by a large set of random initial conditions. 
There is hence no dependence on the initial high-current 
state in our simulation. 

The oscillation regions are shown in Fig. [2j for (a) in- 
creasing and (b) decreasing current and for the angular 
region of interest. The precession frequency varies from 
to about 31 GHz. The critical current for the onset of 
magnetization precession depends strongly on the tilt an- 
gle ft and is reduced almost two orders of magnitude when 
ft increases from almost in-plane (ft = 2°) to perpendic- 
ular (ft = 90°). Both the trend and the quantitative val- 
ues agree well with experiments and other simulations 

@, us mm. 
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FIG. 3: Current-driven hysteresis loops for different ft = 
3° , 8° , 15° , and 20° , showing large hysteresis in both the upper 
threshold current J C 2 and the transition between the static S 
and N states. 



One observes in Fig. [2] that the precession region is 
asymmetric and depends on the direction of the current 
sweep. In the low angle region, both the lower and upper 
threshold currents for precession exhibit hysteresis. We 
define |J c i,+ | and |^c2,+ | as the lower and upper (abso- 
lute) threshold currents for increasing |J|, and similarly 
| «/ci,— | and |J C 2,-| as the corresponding currents for de- 
creasing | J|. As seen in Fig.[2|(a) | J C 2,+ | can be more than 
five times greater than \J C 2,- \ at small ft. The hysteresis 
in J c \ is less obvious in Fig. [21 but will be discussed in 
detail in Fig. [Hat the end of the paper. 

Using Fig. [2j we can now add information about the 
dynamical steady states and their boundaries to the 
phase diagram in Fig.HJ It is noteworthy that the bound- 
aries in the two figures agree. Fig. [2] first confirms our 
assumption that the lack of a stable static state infers 
the existence of a limit cycle in the L region. Secondly, 
it adds a steady dynamic state to the region where our 
eigenvalue analysis only indicated TV; we hence label this 
region N/(L). While Fig. [2] indicates that this state only 
exists for ft > 2°, it does indeed extend all the way to 
the N/(L) — N boundary at about ft = 1°, as was con- 
firmed by choosing initial conditions closer to precession. 
For ft > 2° the S state at P/AP orientation must trans- 
form into L, but below ft = 2° it can simply turn into 
the equivalent TV state at P/AP orientation. Once in the 
N state there is no energetically favorable path to the L 
state. Finally, one realizes that the asymmetry in Fig. [2] 
stems from the selective realization of either the N or L 
state in the N/ (L) region. When approaching this region 
from above in a node state, the systems stays in the node; 
if this region is approached from above in a spiral state, 
the systems enters the dynamic precessional state L. To 
confirm this picture, we simulated minor loops where we 
limited the current sweep to remain in the S/N region 
before reversing the current direction. In this case the 
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a surprisingly rich phase diagram of static and dynamic 
states in zero magnetic field. The coexistence of several of 
these states leads to a number of large hysteretic switch- 
ing behaviors between both static and dynamic states 
and in particular to hysteresis in the threshold currents 
for magnetic precession. 
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FIG. 4: Reduced MR vs. J for ft = 2°. Inset I shows the 
wide-angle orbit of the L state within the S/ (L) region. Inset 
II shows the angular dependence of the hysteretic threshold 
current J c \. 

high current N state is never realized and the S state 
again nucleates a precessional state already at the high 
value of |J C 2,+ | an d not at the much lower |J C 2,-|- 

Not only is the precessional state hysteretic, the two 
static states in the S/N region also exhibit hystere- 
sis. In Fig. [3] we plot the reduced magnetoresistance 
r = (R — Rp)/(Rap — Rp), where Rp and Rap de- 
note the resistance in the parallel and antiparallel con- 
figurations respectively, at four different polarizer angles. 
Close to J = we observe the usual spin torque switch- 
ing region between P and AP states. As the current is in- 
creased we first observe a linear r vs. J region at all angles 
characteristic of the average resistance within the preces- 
sional state L. As precession stops at | J C 2,+ |, r reaches a 
plateau characteristic of the S state located close to the 
north/south poles. At a certain current value (Js^n), 
this state becomes unstable and m switches to its TV 
state close to P/AP alignment. If the current is again 
decreased, m stays within its N state well below Js^n 
and only switches back at a much smaller current, either 
to an S state at Jn^s (ft > 5°) or to an L state at \J C 2,-\ 
(ft < 5°). At about ^=20°, this hysteresis disappears and 
m rotates continuously S <-> N. 

We now finally turn to the hysteresis in the onset cur- 
rent of precession, J c \. In Fig. 0] we plot r vs J in the 
low- J region for ft = 2°. J c \ is clearly hysteretic with 
|J c i,+ |=0.3Jo and | J ci _ |=0.15J . For 0.15 < J/J < 0.3 
the S type P / AP states hence coexist with a precessional 
L state with a very wide cone angle (inset I in Fig. [4]). 
This hysteresis persists at all polarizer angles as shown 
in inset II in Fig. [H 

In conclusion, we have shown by non-linear system 
analysis and magnetodynamical simulations that a spin 
valve with a tilted fixed layer magnetization possesses 
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